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00 ; ABSTRACT 



Context. White dwarfs (WDs) in cataclysmic variables (CVs) are important experimental laboratories where the electron degeneracy 
, , ^ is taking place on a macroscopic scale. Magnetic CVs increase in number especially in the hard X-ray band (> 10 keV) thanks to 

, sensitive hard X-ray missions. 

Aims. From X-ray spectroscopy, we estimate the masses of nearby WDs in moderately-magnetized CVs, or Intermediate Polars (IPs). 
Methods. Using the Suzaku satellite, we aquired wide-band spectra of 17 IPs, covering 3-50 keV. An accretion column model of 
Suleimanov et al. (2005) and an optically-thin thermal emission code were used to construct a spectral emission model of IPs with 
^ , resolved Fe emission lines. By simultaneously fitting the Fe line complex and the hard X-ray continuum of individual spectra, the 

■ shock temperature and the WD mass were determined with a better accuracy than in previous studies. 

' Results. We determined the WD masses of the 17 IPs with statistical fitting errors of < 0.1 - 0.2 Mq in many cases. The WD mass 

I of a recently-found IP, IGR J17195-4100, was also estimated for the first time (1.03^q ,2 ^o)- The average WD mass of the sample 

is 0.88 ± 0.25 Mq. When our results were compared with previous X-ray mass determinations, we found significant deviation in a 
few systems although the reason of this is unclear The iron abundance of the accreting gas was also estimated, and confirmed the 
(N ■ previously reported sub-solar tendency in all sources with better accuracy. 

' Key words, accretion, accretion disks - novae, cataclysmic variables - X-rays: binaries 

(N ■ 
m 

^ ' 1 . Introduction ments of polarization (e.g. lButters et al.ll2009l) . In both types of 

\^ • these binary systems, mass accretion is taking place from non- 

I White dwarfs (WDs) are attractive experimental laboratories degenerate companions, probably low-mass stars, via Roche- 

where the electron degeneracy plays an important role for sup- lobe overflow. Polars do not have an accretion disk. Most IPs 

porting the stars against gravity. An accreting WD is of great im- have accreting disks, while in the inner region, the disk is trun- 

portance in the universe because some of them probably cause cated, and gas almost freely falls onto a WD channeled along its 

y~j Type la supernovae as the WD mass reaches the Chandrasekhar magnetic field. Near the WD surface, a stationary shock stands to 

/\ ' limit, and deliver synthesized heavy elements out to the sur- convert the kinematic energy of the bulk gas motion into thermal 

H rounding space. The mass of a WD is thus one of the critical energy. Because the temperature of the shock-heated gas is typ- 

■ - - 1 physical parameters which determines the fate of the WD. ically > 10 keV, and the density is low, hard X-rays are emitted 

Forming a subgroup of cataclysmic variables (CVs), a mag- 0?^^'^^"^ emission. Then, the heated gas forms 

netic cataclysmic variable (mCV) involves a magnetized WD ^ P^^^ ^^ock region (PSR) with a temperature gradient, wh^ 

as a compact star. mCVs are further classified into polai's and tljgga s descends while i t c2ol^the2^^^ emission (e.g. |Aizu| 

intermediate polars (IPs), according to whether theii magnetic ^ i Frank et al. lEM ^^994; Wu et al. 1994). As a result, 

field strengths are 10^-^ G or 10^-^ G, respectively. Polai-s ai-e so *e total spectrum from flie PSR is obs erved as a sum of multi- 

named after their strong polarization in the optical and infrared temperature emission components (e.g. | lshida et al. || l994|) . 
wavelengths, and sometimes referred to as AM Her type follow - 

ing this prototype (for a review of polars, see ICroppeilll990|) . Although mCVs account for a small fraction (~ 10 - 

The magnetic field of a polar is strong enough for the orbital and 25%) of all CVs, they dominate in the hard X-ray band (> 

spin periods of its WD to be synchronized (Porb - ^'spin)- In IPs, 10 keV), and most (~ 80%) of CVs detected in recent hard 

the two periods are asynchr onous, typically f spin ~ 0.1 Porb (for X-ray surveys with INTEGRAL/ YBIS and Swift are magnetized 

reviews of IPs, see [Patterso n 1994; Hellier 1996), and polariza- (iBarlow et al.ll2006t iLandiet al.|[2 009; Bird et al] |2010l) . includ 



tion is generally weak, although there are some recent measure- ing both known (e.g. iRevnivtsev et al.. ,2004.) and new objects 
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(lAiello et all l2006t lBo"nnet-Bidaud et all l2007t iButters et al.l 



12007': 'Evans et al.ll2008l: ISazonov et al.ll200 8^.'Landi et al."2009') 



A ccording to rece nt X-ray studies (Sazonov et al. 2006; 

IRevnivtsev et al.l l2008). mCVs form an important X-ray source 
population in the lower luminosity range (Lx ~ lO^'-^^^ergs-i), 
while low-mass X-ray binaries with neutron stars and black 
hole binaries dominate in the high-luminosity range. It has also 



Bleach etal. 119721; 


IWorrall et al. 1982; Iwan et al. 
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Kovama etal.1 19861 


kaneda et all ,1997; ,Valinia & Marshall 


1998 


; lEbisawa et al. 


120011 iTanakal 20021 lEbisawa et al. 


2005 


; Kovama et al. 


I2007al; IRevnivtsev & Sazonov! 120071 


Revnivtsev et al. 2009.) and the Galactic Center X-ray emission 


(Kovama et al. 1989ilYamauchi et al.ll 19901 Kovama et al 


199a 


[Muno et al.„2004c 'Kovama et al.l 2007a ; IRevnivtsev et al. 


booi 



I20091 " and references therein) including RS CVn systems, dwarf 
novae, and mCVs, each contributing to different X-ray energy 
regimes. Especially in the hard X-ray band, mCVs are attractive 
candidates because of their spectral resemblance to the GRXE 
(Revnivtsev et al. 2006; Krivonos et al. 2007), including both 
continua and iron-K lines. If we obtain a reliable "averaged 
mCV spectrum" in the hard X-ray band, we can directly com- 
pare it with that of the GRXE. For this purpose, the construction 
of a wide-band spectral model of mCVs and estimation of WD 
masses are both mandatory, because mCV spectra can be well 
expressed as a function of the WD mass as described in Sect. |2] 
Although the mass of a WD is a basic physical quantity, its 
determination is generally not easy, since an uncertainty of the 
inclination angle of the orbital plane degrades the accuracy of 
estimation except for eclipsing systems. However, especially in 
mCVs, X-ray spectroscopy provides an independent method; the 
X-ray spectrum is primarily determined by the shock tempera- 
ture, which in turn reflects the M/R ratio, where M and R are 
the mass and radius of a WD, respectively. Combining the X-ray 
estimate on M/ R with a theoretical WD mass-radius relation, for 
example one bv iNauenbergI (1 1 972l) . 
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we can determine M. 

The X-ray spectroscopic approach is in practice subdivided 
into two methods, the continuum method and the emission 
line method. The former requires a sensitive hard X-ray detec- 
tor since the relevant spectral feature, namely exponential cut- 
off which reflects the shock temperature, appears in the hard 
X-ray range (energy > 10 keV). The latter relies on the re- 
solved emission lines (e.g. from Fe), hence requiring a detec- 
tor with a relatively high energy resolution such as an X-ray 
CCD (charge coupled device) camera. The continuum method 
has been app l ied t o the data taken with Ginga ([ishida 1991; 
Cropper et al.' 1998), RXTE (Ramsay 2000; Suleimanov ( 
2005. .2008). and Swift (Brunschwei ger et al. 2009; Landi 



et all 

2OO9I) iFuiimoto & Ishidal(ll997l) andlEzuka & Ishidal (ll999l) "^ 



plied the emission Une method to the ASCA data, deriving the 

metal abun dances as well. 

Suzaku dMitsuda et al.ll20()7h . which is the fifth Japanese X- 
ray astrophysical observatory, has a wide-band (0.2-600 keV) 
energy coverage with a high s ensitivity, realized by the X- 
ray Imaging Spectrometer (XIS; iKovama et a l. 2007b), which 
uses CCDs and the Hard X-ray Detector (HXD; Takahas hi et alJ 
l2007h which consists of semiconductor detectors and scintillator 
crystals. The hard X-ray response of the HXD allows us to di- 



rectly measure the thermal cutoff of an IP spectrum, and hence 
estimate its shock temperature. At the same time, we can use 
the XIS to resolve emission lines, especially Fe-Ko- lines, into 
three typical components; from almost neutral. He-like, and H- 
like species. Thus, Suzaku is expected to provide significantly 
advanced information on the masses of mCVs. 

In the present paper, we report masses and iron abundances 
of IPs estimated by analyzing wide-band Suzaku spectra of 17 
objects. Compared to IPs, hard X-ray emission from polars 
is weaker because cyclotron emission induced by the stronger 
magnetic field reduces their plasma temperatures. It is also dif- 
ficult to incorporate the effect of strong magnetic field into a 
spectral model. Therefore, we only analyzed IPs. 

In Sect. |2] we describe a numerically calculated emission 
model which is used to derive the mass estimations. Target infor- 
mation and data reduction procedure are given in Sect. [3] Section 
|4]deals with the WD mass estimation. In Sect.jSj we compare our 
results with previous reports. 

2. Modeling of X-ray emission from an accretion 
column 

2.1. Previous works 

A number of works have been performed to model X-ray spec- 
tra emitted from PSRs of accretion columns in mCVs, and to 
estimate the PSR temperatures by comparing the model spec- 
trum with observations. The temperatures were then translated 
into the M/R ratio of their WDs. As mentioned in Sect. [Tj these 
works can be divided into continuum and emission line methods. 

In the continuum method, the vertical structure of the ac- 
cretion column (e.g., profiles of bulk velocity, temperature, and 
density) is calculated based on plausible physical assumptions, 
following which an X-ray spectral model can be constructed 
by convolving thin thermal emission models with the plasma 
profile. The M/R ratio is a free parameter of the synthesized 
model spectrum, and should be det ermined thr ough its com- 
parison with an observed spectrum. Aizul (119731) was the very 
first who explored this theoretical method. By neglecting the 
gravitational a ttraction inside accretion columns. IWul(ll994l) and 
IWu et al.l (1199 4) presented profiles of plasma parameters in one- 
dimensional plane-parallel PSR in closed-integral formulae by 
considering plasma cooling via bremsstrahlung and cyclotron 
em ission. 

Ilshidal (Il99lh analyzed Ginga data of several IPs using an 
isothermal assumption to a PSR, and reported WD mass estima- 
tions (or lower limits) for the first time by fitting spectra with 
an isothermal emission model. By comparing the isothermal 
emission model with their more realistic spectral calculations, 
IWu et al.l (I1995 I) gave correction factors to th e W D masses esti- 
mated bv llshidal(1199ll) . ICrop per et al.'('l998') and'Ramsa v et all 
( I99I) also used the model of Wu ( 1994) and fitted data obtained 
with Ginga, ASCA, and RXTE. 

To c onstruct a more p hysic ally realistic modeklCropper et al.l 
(119991) . iRamsavl (120001) and ISuleimanov eTaP (120051) calcu- 



lated the accretion column structure by taking into account 
the gravitational attraction, and found that the inclusion of the 
gravity slightly reduces the WD mass estimates espe cially in 
systems with higher WD masses (e.g. > 1.0 Mq; iRamsavl 
2000). Brunschweiger ejal] (|2009|) also applied the model of 
ISuleimanov et al.i (2005) to the data of 22 IP systems ob- 
served with the Swift/BAI, and gave mass estimates. Some au- 
thor s considered more sophis ticated emission models; for exam- 
ple, ISuleimanov et al.l (l2008l) considered the effect of Compton 
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up-scattering of photons by post-shock plasmas although it 
was found to change the WD mass estimates only slightly. 
[Canalle et al. (2005) calculated PSR temperature and density 
profiles by taking into account the dipole magnetic field struc- 
ture, and showed that this increases estimated WD masses 
mainly in high mass (e.g. ~ 1 Mq) systems. 

As the fi rst application of the emission l i ne m ethod, 
llshida & Filiim oto (1995) and iFuiimoto & Ishidal (11997 ) ana- 
lyzed the ^5CA/SIS data of EX Hya. They thus estimated 
the plasma temperature behind the shock by comparing ob- 
served line-intensity ratios with those theoretically predicted. 
Employing the PSR structure of Aizu (1973), they obtained a 
WD mass estimate of O.48^y ygM0, which agreed with that of 
0. 49 + 0.03Mo d erived with optical spectroscopic observations 
by|Hellier'('l996'). Using_the sa me line -ratio model, Helli er et al.l 
([1996), Fujimoto (1991), and lEzuk a & Ishida (1999) investi- 
gated WD masses in other systems, giving comparable results 
with those obtained with the X-ray continuum methods. 

2.2. The constructed model 

To make the best use of the capability of Suzaku (Sect. 1), we 
constructed a spectrum model by improving the PSR model of 
[Suleimanov et al. (2005) in two ways. One is to include the Fe-K 
emission lines, while the other is to improve the plasma cooling 
function, which is used when solving hydrodynamical equations 
in the high-temperature regime (> 10**-^ K; see the following 
section). Below we describe how we constructed our model. 

As illustrated in Fig.[T] a one-dimensional cylindrical accre- 
tion column, or PSR, is assumed. Channeled by the magnetic 
field, gas is considered to fall into the accretion column with its 
supersonic free-fall velocity. At a certain height, a strong shock 
converts some of the bulk motion of the gas into thermal en- 
ergy in the PSR. While the shock-heated gas falls subsonically 
down along the magnetic field lines, it cools by emitting X-rays. 
Finally, the gas is assumed to softly land onto the atmosphere of 
the WD. Because the magnetic fields of IPs are not as strong as 
those of polars, the dominant cooling process of the heated gas is 
optically thin thermal X-ray emission (free-free plus free-bound 
plus bound-bound), and t he cyclotron coo ling can be considered 
negligible (e.g. Eq.(2) of IWu et aP 1 1 995h . Including these, the 
assumptions made in the present model are summarized as 

- 1-dimensional cylindrical accretion 
(curvature of magnetic field neglected) 

- temperature before the shock to be zero 

- gas to be ideal 

- a strong shock 

- cooling via free-free, free-bound, and bound-bound process 
(cyclotron cooling to be negligible) 

- gravitational force in the PSR taken into account 

- soft landing of cooled gas onto the WD surface 

These are the same assumptions as those used in 
ISuleimanov et al.l (2005) and the negligible magnetic field 
case of lCropper et aL (il999i) . 



2.3. Calculation of the model 

With these assumptions, the accreting plasma flow can be de- 
scribed with hydrodynamical equations; conservations of the 
mass, the momentum, and the energy. Since these equations 
and initial conditions assumed in the numeric al calculation have 
been already published by several authors (e.g. lSuleimanov et al.l 



120051: ICropper et al.ll 1 999h . we do not recall t hem here, and fol- 
low the same formalization and notation as ISuleimanov et al.l 
(2005). Below, we describe major changes and some attentions 

related to the model calculation; 

In the model used in ISuleimanov et al.l (l2005h and 
IBrunschweiger et al] ( 2009h. the p l asma c ooling function A was 
taken from 'Suthe rland & Dopital (119931) . Although we could 
integrate the optically-thin emission code (apec; see below) 
to construct a cooling fun ction, we also use the result of 
[Sutherland & Dopital (Il993h here. This is because spectra cal- 
culated with the emission code only cover £ < 50 keV cur- 
rently, and this would artificially reduce cooling rates in higher 
plasma temperatures (e.g. IcT > 10 keV). However, for rather 
massive WDs, we need to spe cially treat a higher temper ature 
regime which is not covered bv lSutherland & Dopital (Il993 l) (i.e. 
T > 10^^ K). In this regime, we extended the cooUng function 
using a thermally averaged relativistic free-free Gaunt factor tab- 
ulated by Nakagawa et al., (1987) . We averaged the Gaunt factor 
over the frequency separately for hydrogen and helium, and then 
converted the averaged values into the cooling rate (or energy 
loss rate) in terms of the total po wer per unit volume emitted by 
thermal bremsstrahlung (see e.g. lRvbicki & Lightmanl[T979l) . 



dW 
dVdt 



1.426 X lO^^^neWiZ^gB erg cm"^"* s"'. 



(2) 



where We, "i, Z, and are electron and ion number densities, 
charge of ion, and the averaged Gaunt factor, respectively. Then, 
by rescaling the obtained cooling rate so that it matches that of 
Sutherland & Dopita ( 1993) in T < 10**^ K, we constructed the 
cooling rate up to T < 10^ " K. The extrapolated cooling function 
exceeds a simple linear extr apolation of the cooling function by 
[Sutherland & Dopital (1 19931) bv a factor of < 2 . 

When a fraction / of the PSR cross section to the WD sur- 
face area is set, the specific accretion rate a can be related to 
the total accretion rate M via M = AnR^^af. Since the lumi- 
nosity of a PSR is determined by liberated gravitational energy 
(oc M), the product af can be regarded as the model normaliza- 
tion. Although a and / can be free parameters in general (see e.g. 
Cropper et al. 1999; Suleimanov et al. 2005), and the numerical 
solution changes as a varies, we fixed them at a = 1 .0 g cm"^ s"' 
and / - 0.001 (0.1% of the WD surface area) for the sake of an 
easy comparison with the previous reports. To rescale the model 
spectrum, a normalization parameter, which linearly adjusts the 
normalization of the model calculated with the fixed a and /, was 
introduced. Eff'ects of changing a on our result are discussed in 
Sect. l5A2l 

After setting A and a, finding the solution of the three con- 
servation equations is reduced to solving an initial value problem 
of ordinary different ial e quations. As explai i ied fo r example in 
ICropper etaP (119991) and ISuleimanov et al ] (120051) . this can be 
practiced with the so-called shooting algorithm. Each trial starts 
from an assumed shock height z - Zq, where the initial param- 
eters are set according t o the strong shock l imit of the Rankine- 
Hugoniot relation (Suleimanov et al.ll2005b . The lower limit of 
the calculation, z - Rwd, is specified by the WD mass-radius re- 
lation of Eq.([T]i- At the WD surface, another boundary condition 
of the bulk gas speed v = is imposed to fulfill the soft-landing 
assumption. In the algorithm, we first solve the differential equa- 
tions from a trial initial value of zo (for example, 1.01 Rwd), and 
see how the final state of each quantity at z = /?wd deviates from 
the conditions assumed above. By taking into account the devi- 
ation as a feedback to a new initial zo, we solve the equations 
again. Until the feedback becomes small enough, below ~ 10"^ 
of derived zq, this procedure is repeated. Practically, this method 
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converged in several to a few tens of iterations at each value of 

MwD- 

We executed the calculation over a range of 0.20 Mq < 
MwD < 1-39 Mq with a step of 0.01 Mq. Figure |2] shows the 
dependency of the shock height and the shock temperature on 
the WD mass by interpolating data points with spline functions. 
In the higher WD masses, the gas velocity at the shock height 
becomes higher, and therefore the gas density falls for a given 
accretion rate. This causes a decrease in the plasma cooling rate 
an d cooling time, yielding taller shock heights. 

ICropperetaTI ( 119991) showed that an inclusion of the grav- 
itational attraction changes WD mass estimations especially in 
high mass (~ 1.0 Mq) systems. This can be confirmed with the 
comparison of the profiles of shock temperature calculated with 
and without gravitational attraction shown in Fig. |2] For exam- 
ple, at MwD -1-2 Mq, no-gravity profile overestimates a shock 
temperature and following Mwp by ~ 10% co mpared with the 
present result. Note that in ICropper et al.l (11999) . they used an 
analytic free-free cooling rate ( Wull 19941) ? which is described in 
their Appendix A. Therefore, the resulting shock temperatures 
(their Fig. 1) are slightly different from our values (Fig.|2]i. 

Based on the calculation results, we convolved profiles of 
the temperature and the density along the accretion column 
with (single-temperature) spectra calculated by the optically- 
thin emission code apec (ISmith et al.l 1200 ll) to obtain a total 
multi-temperature emission spectrum. Since we employed apec 
spectra instead of the si mple bremsstrahlung continuum which 
ISuleimanov et al.l (12005 ) utilized, our model includes atomic 
(particularly Fe) emission lines, which can be an important sig- 
nature of the plasma temperature in the present spectral analysis. 
In this way, we prepared a series of model spectra for each Mwd 
by changing the metal abundance of apec from 0. 1 Z© to 2.0 Zq 
with a step of 0.1 Zq. 

Because the currently available version of apec does not cal- 
culate the spectrum above £ > 50 keV, we can miss the expo- 
nential cutoff feature, which is crucial in the calculation of high 
temperature layers (kT > 20 - 30 keV). Although this energy 
range is not directly employed in our data analysis, the model 
shapes in that range are considered important, because the HXD 
has the broad spectral redistribution function, and photons with 
£■ > 50 keV sometimes produce signals in significantly lower 
energies. Therefore, above £ > 50 keV, we substituted spec- 
tra estimated by meka model (which still provides a spectrum in 
that energy range), with its normalization set 5.5% smaller than 
that of the corresponding apec model. This renormalization is 
done to make the two spectra match smoothly at £ = 50 keV 
under high temperatures where this substitution is of particular 
importance. In this manner, the model was calculated over the 
0.1 - 100 keV band. Note however that we only used it only 
above 3 keV in spectral fitting to avoid modeling of complex ab- 
sorption in lower energies, which is less important for WD mass 
estimations (see Sect. 14. lb . 

Near the WD surface, where the density of a PSR plasma 
rapidly increases, the optically-thin assumption will be invalid, 
and optically-thick emission will emerge. Although this can alter 
a PSR spectrum in the lower energy range (e.g. £ < 0.1 keV), we 
simply used apec because our analysis (Sect.O utilizes spectra 
in higher energies {E > 3 keV). 

Finally, a set of total spectra calculated in this way were con- 
verted to a so-calle d "local mode l" file of the spectral analy- 
sis software Xspec (lArnaudlll996h . The free parameters of the 
model are the WD mass, the metal abundance, and the model 
normalization. In this local model, spectra for masses and abun- 
dances that are in between the calculation grid points (see above) 
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Fig. 1. Assumed geometry of the accretion column of an IP. 
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Fig. 2. Results of the numerical solutions for the shock height 
from the WD surface (thick solid line) and the shock temperature 
(thin solid line), shown against the WD mass. For comparison, 
the dashed line shows the shock temperature calculated by as- 
suming no-gravity in the PSR. 



are generated by interpolating spectra of the surrounding four 
grid points. 

In Fig. [3] we present normalized profiles of the temperature 
and density for an IP with a WD mass of 0.7 Mq, for comparison 
with the results of iSuleimanov et al.i (|2005). The shock height of 
0. 013 RwD, which we deriv ed, is somewhat different from that 
of iSuleimanov etal.l (l2005h . 0.018 R^^u■ This difference might 
arise from difference of implementation of the algorithm. In the 
present study, we do not investigate this difference since the 
shock temperatures in the two cases differ only by ~ 0.5%. 



3. Observations and data reduction 

We assembled our IP sample of 17 systems observed with 
Suzaku. Table [T] lists information of their Suzaku observations 
together with their background-subtracted count rates from each 
detector. 
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Fig. 3. Temperature and density profiles of a PSR calculated 
for an IP of Mwp - 0.7 Mq. For easy comparison with 
[Suleimanov et al.l (l2005h . the plot employs the same parameter 
values and the same style as their Fig. 2. Labels give the WD ra- 
dius calculated using the mass-radius relation of Eq.([T]l, as well 
as the shock height (zq-Rwd) and the shock temperature (fermax) 
from Fig.|2] 



The XIS consists of four X-ray imaging charge coupled 
device (CCD) cameras as a focal plane detector of the X-ray 
Telescope (XRT; Serlemitsos et al. 2007). Three (XISO, 2, and 
3) use front-side illuminated (FI) CCD chips, while the other is 
back-side illumin ated (BI; XISl). The e nergy coverage is nomi- 
nally 0.2-12 keV (iKovama et al.ll2007bh . Since 2006 November, 
one of the FI sensors, XIS2, has not been operational, and there- 
fore, we utilize data taken with the other three XIS sensors for 
observations performed after that date. 

The HXD is a non-imaging collimated X-ray detector that 
consists of two main detection parts; silicon /9-intrinsic-n diodes 
(hereafter PIN) and gadolinium silicate (GSO) crystals. They 
cover the energy range of 10 - 70 keV (PIN) and 40 - 600 keV. 
In the present analysis, we only utilize PIN data because signals 
from IPs were hardly detected by GSO. 

Since the sample has been constructed by combining differ- 
ent observations, it is heterogeneous and incomplete, although 
it can be regarded as a flux-limited sample in the hard X- 
ray band. To perform a detailed spectral analysis, we only se- 
lected targets with enough statistics in the hard X-ray band; 
sources with background-corrected HXD count rates higher than 
0.04 counts s"' are used in the present analysis. This rate is 
equivalent to 10 % o f the typical HXD bac kground count rate 
dKokubun et al.ll2007l: iFukazawa et al.ll2009l) . or ~ 1 mCrab in 
terms of flux. With these criteria, AE Aqr, GK Per, and the IP- 
like object SAX J1748.2-2808 (Nobukawa et al. 2002) were ex- 
cluded from the sample, although the XIS clearly detects their 
spectra. 

The typical energy bands utilized in the present study are 
3 - 12 keV and 12 - 50 keV for the XIS and the HXD (Table 
[U. We limited ourselves to use the XIS data above 3 keV since 
lower energy spectral shape are afifected by the intrinsic multi- 
column absorption; the time-averaged spectrum is integrated 




10 

Energy (keV) 




1 10 

Energy (keV) 

Fig. 4. A calculated total spectrum from a PSR of an IP 
with a WD mass of 0.7 Mq and a gas abundance of 1.0 Zq, 
plotted in units of (a) photons cm~^ s~' keV"' and (b) 
keV^ cm"^ s"' keV"'. The calculation assumed a distance of 
100 pc to the WD. 



over multiple azimuthal angles, and thus including spectra ab- 
sorbed by different column density. The higher floor bound of 
the XIS in some sources are explained in the section below. 
In some observations, especially in that of FO Aqr, we applied 
higher floor bound for the HXD data since an enhanced electri- 
cal noise contaminated the data therein. The HXD data above 
40 keV are discarded when source signals were below 5 % of 
the background. 

All observations were performed in the normal operation 
mode of the XIS and HXD. Obtained data were processed with 
the Suzaku off-line processing software version 2. We utilized 
the analysis software package HEASOFT of version 6.6.2 to- 
gether with a detector calibration database (CALDB) version 
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2009-05-11. We excluded the XIS and HXD/PIN data taken 
when the target is occulted by the Earth, or the spacecraft was 
in the South Atlantic Anomaly. To reduce the contamination 
from scattered Solar X-rays, we also removed the XIS data taken 
while the target elevation angle above the Sun-lit Earth rim was 
less than 20°. We applied the same event filtering criteria as 
used for creating nominal cleaned event files, as described at the 
Suzaku websitqH 

In the XIS data analysis, we extracted time-averaged spectra 
of the targets using circular regions with radii of 2'. Background 
spectra were also extracted from the same observation data, ex- 
cluding apparent point sources if any. This background spec- 
trum contains both non-X-ray background (NXB) and diffuse 
X-ray background (DXB). Because the signal-to-noise ratio is 
sufficiently high for all targets, the background subtraction ig- 
nored the vignetting efi'ect of the XRT (Serlemitsos et al. 2007). 
We generated redistribution matrices and auxiliary response files 
of the XI S/XRT using xisrmf gen version 2009-02-28 and xissi- 
marfgen (llshisaki et alJl2007l) version 2009-01-08, respectively 
(both contained in HEASOFT). In the spectral analysis, we 
summed the spectra taken by the Fl sensors, and used a detector 
response averaged over them. 

We accumulated a time-averaged PIN spectrum of each tar- 
get, and subtracted the non X-ray background events extracted 
from the NXB event file, which in turn was produced from 
a modeling method by Fukazawa et al. (2009) and distributed 
via the PIN NXB release websit^. The version keyword of the 
NXB file we used is METHOD =LCF1TDT, i.e. so-called tuned 
ba ckground. We subtracted the DXB component simulated af- 
ter Boldt (1987), and subtracted it from the observed spectrum. 
This procedure is summarized in the website of NASA Goddard 
Space Flight Centefl Since the present sample objects are lo- 
cated off the Galactic plane, we neglected contributions from the 
GRXE. When fitting a PIN spectrum, we employed an appro- 
priate redistribution matrix file contained in CALDB, depending 
on the pointing mode and the observation date. Because the ob- 
servation of XY Ari was performed at an off-axis position of the 
HXD, the effective area of PIN decreased by ~ 20%; we took 
this effect into account by applying an additional auxiliary file 
calculated with hxdarfgen. 

In some observations, XIS images contained weak point 
sources other than target IPs. These sources could contaminate 
the data of PIN because PIN has a wider field of view ( 1° x 1° bot- 
tom to bottom). However, their contributions in the PIN energy 
band are estimated to be < 2 - 3% in XY Ari, and much lower in 
other observations. Therefore, we simply neglected them in the 
present analysis. 



4. Spectral analysis and results 

4. 1 . Fitting to tlie XIS/HXD spectra witii tlie partial covering 
modei 

Using the tabulated model (Sect. |2l hereafter ip_psr), we 
performed chi-square fitting to the spectra of individual IPs. 
The lower energy band (< 3 keV) was ignored because anal- 
yses in this band are hampered by multi-col umn and ion- 
ized absorption caused by pre-shock gas (e.g., ICropper et alJ 
[T999tlde Martino et al.ll2001 ). To account for the intrinsic multi- 
column absorption, we utilized a model with a partial-covering 



' http://www.astro.isas.jaxa.jp/suzaku/process/v2changes/criteria_xis 

^ http://www.astro.isas.ac.jp/suzaku/analysis/hxd/pinnxb/ 

^ http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/pin_cxb.html 



(PC) absorber expressed by vphabs (available in Xspec). A 
Gaussian was added to the model to express the neutral (or low 
ionized) Fe Ka emission line at ~ 6.4 keV. In addition, the entire 
model spectrum was subjected to interstellar absorption using 
phabs. The total model expression is thus 

constant x phabs x 
(Cpc ■ vphabs ■ ip_psr 
+(1 - Cpc) ■ ip_psr + gaussian), 

where constant is a renormalization factor for the two detec- 
tors (see Suzaku Memo 2008-0^3); fixed at 1 for the XIS, and 
1.16 (XlS-nominal attitude) or 1.18 (HXD-nominal attitude) for 
PIN. Cpc denotes a partial covering fraction. The two ip_psr 
components are constrained to have a common set of parameters. 
The Fe abundance of vphabs was equalized to that of ip_psr, 
considering that the X-ray absorbing and emitting gas must have 
the same chemical composition. 

First, we fitted the wide-band spectra using all XIS data start- 
ing from 3 keV. The HXD spectra were fitted over the energy 
ranges listed in Table [1] The table also gives the shock tem- 
peratures, which are related to Mwd through Fig. |2l Hereafter, 
the quoted errors are at 90% confidence levels unless otherwise 
noted. Note that the WD mass of IGR J17195-4100 was esti- 
mated for the first time as far as we know. 

In the present fitting, we considered fits with a null- 
hypothesis probability higher than 1% as acceptable. Based on 
this criterion, acceptable fits were obtained with this partial cov- 
ering model in 1 1 systems; FO Aqr, EX Hya, NY Lup, V2400 
Oph, V1223 Sgr, and RX J2133 are remainders. The best-fit 
parameters are listed in Tableland Figs. |5] and I A. II show the 
spectra and the best-fit models of individual IPs. For the remain- 
ing six objects, we listed results obtained after prescriptions de- 
scribed in the paragraphs below. The average WD mass and the 
mean iron abundance associated with their Icr standard devia- 
tions were 0.88 + 0.25 Mq and 0.38 + 0.14 Zq, respectively. 

When fitted with the partial covering model, EX Hya showed 
large bipolar residuals in the energy b and slightly below t he H- 
like Fe Ka line. As reported by Hel lier & Mukail (|2004 with 
ChandrafHETG data, the highly ionized Fe lines from EX Hya 
have complex shapes with considerable contribution from unre- 
solved dielectric satellite lines. Although this might be a cause of 
the residual, we leave this speculation as a future task because it 
is difficult to diagnose the line structure with the energy resolu- 
tion of the XIS. In the present study, we therefore put the original 
fitting result obtained with the partial covering model in Table|2] 
The estimated WD mass did not change within errors even if the 
H-like Fe emission line was ignored in the fitting since the ther- 
mal cutoff is clearly detected in the PIN spectrum, and it strongly 
constrains the shock temperature (and hence a WD mass). 

In NY Lup, V1223 Sgr, and RX J2133, positive residuals 
were observed near the lower bound of the XIS data. We sus- 
pect that these residuals reflect the complex absorption structure 
in the lower energy spectra, which are not as so important when 
determining shock temperatures and WD masses. Therefore, we 
gradually changed the lower bound of the XIS spectrum, and 
found that fits become acceptable when the lower bound ener- 
gies listed in TablelU 3.5 keV for NY Lup, and 4 keV for V1223 
Sgr and RXS J2133, were employed. These changes did not sig- 
nificantly affect the best-fit parameters, and therefore we simply 
regarded the resultant parameters of these fittings with the nar- 
rower energy bands as the current best estimations, 
html 

ftp://legacy.gsfc.nasa.gov/suzaku/doc/xrt/suzakumemo-2008- 
06.pdf 
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System 


Coordinate" 


Start Time 


Exp." 




Enerj 


;y Band*^ 




Rate" 


Obs. ID 




RA 


Dec 






ks 


XIS (keV) 


HXD (keV) 


XIS 


HXD 




FO Aquarii 


334.481 


-8.351 


2009 


06-05 08:14 


33.4 


3.0- 


12.0 


17.0 - 


40.0 


0.95 


0.08 


404032010 


XY Arietis 


44.036 


19.442 


2006 


02-03 23:02 


93.6 


3.0- 


12.0 


12.0- 


50.0 


0.65 


0.07 


500015010 


MU Camelopadahs 


96.318 


73.578 


2008 


04-14 00:55 


50.1 


3.0- 


12.0 


12.0 - 


40.0 


0.27 


0.06 


403004010 


BG Canis Minoris 


112.871 


9.94 


2009 


04-11 12:11 


45.0 


3.0- 


12.0 


12.0- 


40.0 


0.51 


0.10 


404029010 


V709 Cassiopeiae 


7.204 


59.289 


2008 


06-20 10:24 


33.3 


3.0- 


12.0 


12.0- 


50.0 


1.10 


0.16 


403025010 


TV Columbae 


82.356 


-32.818 


2008 


04-17 18:00 


30.1 


3.0- 


12.0 


12.0- 


50.0 


1.21 


0.20 


403023010 


TX Columbae 


85.834 


-41.032 


2009 


05-12 16:19 


51.1 


3.0- 


12.0 


12.0- 


40.0 


0.29 


0.04 


404031010 


YY Draconis 


175.896 


71.703 


2008 


06-15 18:37 


27.4 


3.0- 


12.0 


12.0- 


40.0 


0.85 


0.10 


403022010 


PQ Geminorum 


117.822 


14.74 


2009 


04-12 13:46 


43.2 


3.0- 


12.0 


13.0- 


50.0 


0.56 


0.09 


404030010 


EX Hydrae 


193.107 


-29.249 


2007 


07-18 21:23 


91.0 


3.0- 


12.0 


12.0- 


33.0 


2.12 


0.14 


402001010 


NY Lupi 


237.061 


-45.478 


2007 


02-01 15:17 


86.8 


3.5- 


12.0 


12.0- 


50.0 


0.79 


0.19 


401037010 


V2400 Ophiuchi 


258.146 


-24.247 


2009 


02-27 11:42 


110 


3.0- 


12.0 


12.0- 


50.0 


2.60 


0.21 


403021010 


AO Piscium 


343.825 


-3.178 


2009 


06-22 11:50 


35.6 


3.0- 


12.0 


12.0- 


50.0 


1.21 


0.15 


404033010 


V1223 Saggitarii 


283.759 


-31.163 


2007 


04-13 11:31 


46.2 


4.0- 


12.0 


13.0- 


50.0 


2.38 


0.52 


402002010 


RX J2133.7+5107 


323.432 


51.124 


2006 


04-29 06:50 


62.8 


4.0- 


12.0 


12.0- 


50.0 


0.83 


0.15 


401038010 


IGR J17303-0601 


262.59 


-5.993 


2009 


02-16 10:09 


27.7 


3.0- 


12.0 


13.0- 


50.0 


0.56 


0.16 


403026010 


IGR J17195-4100 


259.898 


-41.015 


2009 


02-18 11:03 


26.9 


3.0- 


12.0 


12.0- 


50.0 


1.02 


0.15 


403028010 



" In units of degree (J2000 equinox). " Net exposure. Energy band used in the spectral analysis (Q. Background-corrected count rates in 
units of counts s"' integrated over the noted energy bands. 



In FO Aqr and V2400 Oph, the partial covering model did 
not give an acceptable fit even when the employed energy ranges 
were changed. Instead, the fits improved significantly to an ac- 
ceptable level when an additional absorbed ip_psr component 
was incorporat ed. Since these sources a re supposed to have high 
accretion rates (ISuleimanov et al.ll2005h . we consider that taking 
this additional absorbed component into account is reasonable. 
Therefore, we estimated WD masses in these sources using this 
modified model. The column densities of the additional absorbed 
component are also listed below those of the primary PC model 
in Table 12] 

In systems with special accreting and viewing geometries, 
the Fe emission lines (6.7 keV) could be artificially enhanced 
or suppressed by "r esonance trapping" of the Fe-K line photons 
(iTerada et al.ll200l\ because the optical depth of the process is 
higher than 1 (e.g. 25 by Done et al. 1995). This could affect our 
abundance measurements, and in turn, those of Mwd- However, 
this effect is only of a importance when the column is observed 
either pole-on (the Fe-K line enhanced) o r side-on (suppress ed) 
through a large fraction of the spin phase ( Terada 



a (suppress es 
etal.ll2001h . 



Because our 17 IPs are more or less pulsating at their spin 
periods, they are not likely to be either pole-on or side-on, and 
the enhancement or suppression, if any, can be averaged out by 
the WD spin (changing viewing angle). We therefore consider 
that this effect is negligible for the present study, and that the 
derived sub-solar abundances reflect a true property of the ac- 
creting gases. 

4.2. HXD/PIN-only fitting 

We also fitted the HXD/PIN spectra alone to confirm the result 
of the partial covering model at the energy band where even the 
strongest absorption little changes the spectral model. In the fit- 
ting, we fixed Cpc - 0. 

The results of this analysis are compared with those from 
the wide-band fitting (Sect. 14.11 ) in Fig. |6] Large y-axis errors 
are simply due to lower statistics in the PIN data. The masses 
derived with the two methods thus agree with each other within 
errors. In particular, we do not find any systematic differences 



1.4 



1.2 



0.8 



0.6 



0.4 
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IGRJl 
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BG 

N¥-'' 



7303 



PCI V709 



XY 



TX 



0.2 0.4 0.6 0.8 1 1.2 
M\yQ (This work, Partial Covering) 



1.4 



Fig. 6. Comparison of the WD masses estimated from the par- 
tial covering model using the XIS and HXD/PIN data, with those 
determined by the HXD/PIN data alone ignoring the partial cov- 
ering. Source names are abbreviated. 



between the two sets of estimates. This self-calibration ensures 
the consistency of the continuum and line methods. 

5. Discussion 

5. 1 . Comparison witli previous reports : WD masses 
5.1 .1 . Based on X-ray data 

Table|3]compares the present Suzaku XIS -H HXD results with the 
WD masses reported by the previous studies in the X-ray band. 
In Fig.|2]we plot correlations between our WD masses and those 
from the RXTE /PC A+HE XTE da ta by Suleimanov et al. (2005j) 
and the Swift/BAI databv lBrunschweiger et aL(i2009.) : the spec- 
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FO Aquarii FO Aquarii (deconvolved) 




Energy (keV) Energy (keV) 



XY Arietis XY Arietis (deconvolved) 




Energy (keV) Energy (keV) 



Fig. 5. Observed XIS and HXD spectra (black crosses) of individual IPs compared with the best-fit model (solid line) . Left panels: 
Raw spectra, wherein the detector responses are not removed. The lower panels show the data-to-model ratio. Right panels: The same 
spectra, shown after removing of the detector responses and multiplied by E-. Equivalent to the vFy presentation. The discontinuity 
seen in the model curve at 12 keV is due to the renormalization factor, which was introduced to adjust relative effective areas of the 
XIS and PIN (see text). Remaining figures are shown in Fig. A.l of Appendix A. 



tral models and the energy band used in these three studies are 
very similar to one another. 

In Figs.|2](a) and (b), our results derive higher WD masses, 
especially for massive (> 1 Mq) systems, although rough cor- 
relations are seen among the three measurements. For a few 
sources, for example V709 Cas, BG CMi, and PQ Gem, this 
discrepancy amounts to ~ 0.3 - 0.4 Mq. Reasons for these dif- 
ferences are unclear at pr esen t. Unfortunately, the W D masses of 
ISuleimanov et al.l (l2005h and lBrunschweiger et al.l ( [2009 ) show 
the same level of disagreement as seen in these figures, although 
the two studies use the same spectral model calculated in the for- 
mer This might indicate that the results suffer from rather large 
systematic uncertainties, which could exceed the quoted statisti- 
cal fitting errors. We are aware that our analysis is not free either 
of possible systematic uncertainties (Sect. I5.4l i involved in the 
detectors and the spectral model. Yet our result may be less sub- 
ject to various systematic effects, because we jointly utilized the 
Fe-K line spectroscopy and the hard-band continuum shape, and 
the wide-band fitting result is consistent with that from the PIN- 
only fitting (Fig.|6]l. 

From X-ra y light curve analysis in the eclipsing IP XY Ari, 
iHellied (Il997h estimated Mwd to be 0.91 - 1.17 Mq. Our re- 
sult, MwD = 1-14^0 [4 Mq, nicely agrees with this value, and 



this agreement provid es a convincing calibration to the present 
methodology (see also lBrunschweiger et al.ll2b09l) . 

5.1.2. Other methods 

It is also meaningful to compare the derived WD masses with 
those estimated from binary kinematics based mainly on op- 
tical/IR spectroscopy. Using references in the CV catalog by 
Ritter & Kolb (2003), we collected reports on Mwd of IPs, and 
then plotted them as shown in Fig. |8] Note that most of the re- 
ports employed an assumption of a Roche-lobe filling Zero Age 
Main Sequence star to estimate radius and mass of the secondary 
(non-degenerate) stars, and this may be a cause of systematic un- 
certainties in their WD mass estimations. 

All masses agree (with rather large errors) except for EX 
Hya. For this object, WD masses determined spectroscopically 
in optical/IR wavelength have been somewhat controversial. The 
difference of the WD masses are mainly caused by different ra- 
dial velocity estimation s. Based on Nal and Call line profiles, 
Beu ermann & ReinschI (2008) estimated the amplitude of a ra- 
dial velocity of the secondary (non-degenerate) star K2 to be 
432.4 + 4.8 km s~' by assuming that of the primary WD K \ 
of 59 + 3 km s ' dBelle et al.1 120031: iHoogerwerf et all |2004 . 
From these values, they obtained Mwd = 0.790 + 0.026 Mq. On 



Table 2. 

levels. 
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The best-fit spectral parameters obtained with the partial-covering absorber model. The errors are at 90% confidence 
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1.11(457) 


1.59 


2.73 


1.02(291) 


1.03 


1.79 


0.98(350) 


2.24 


3.18 


0.99(382) 


4.12 


5.42 


1.08(383) 


5.08 


5.92 


1.18(341) 


1.23 


1.18 


1.07(340) 


3.77 


2.74 


1.17(386) 


2.49 


3.11 


1.22(358) 


9.10 


3.56 


1.19(307) 


3.29 


6.03 


1.12(457) 


5.00 


5.99 



1.03(595) 


4, 


.96 


4.06 


1.19(333) 


9, 


.52 
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2, 


.10 
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2, 


.07 
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3, 
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"The absorption column density of fully-covering absorbent, i' The absorption column density of the partial-covering absorbent (see text). 



"^The partial-covering fraction. 'i2 - 10 keV flux in units of 10 
density of the additional absorbed component (see text). 



erg cm 2 s 1. ''12 - 40 keV flux in units of 10 " erg cm 2 s 1. ^The column 



Table 3. Estimated WD masses compared with previous reports in the X-ray band. The mass ranges associated are all 90% confi- 
dence intervals. 
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MU Cam 
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0.85 ±0.12 
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1.09(> 0.94) 




V709 Cas 


J 22+0.05 


0.96 ± 0.05 
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TV Col 
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TXCol 


67+^-^^ 


0.67 ±0.10 


0.70 ± 0.30 


74+°-"^ 

"■'^-0.05 


0.48!SJ;S^ 


YY Dra 


0.50 ±0.11 


0.75 ± 0.05 








PQ Gem 




0.65 ± 0.09 


0.65 ± 0.20 




1.29(> 1.17) 




EXHya 


0.42!|| 


0.66 ±0.17 


0.50 ± 0.05 


44+°-"3 

"■^^-0.03 


46+" °'' 

"■^"-0.04 


o.48!S:i° 


NY Lup 


1 15+"-o^ 


1.09 ± 0.07 










V2400 Oph 


62+0-06 
-fl.QS 


0.81 ±0.10 


0.59 ± 0.05 


0.71+°-"^ 




0.68!»| 


AOPsc 


61+"* 

75+"-°^ 
-fl.95 


0.55 ± 0.06 


0.65 ± 0.05 


60+^-^^^ 

"■""-0.03 


56+" 16 
"■-'"-0.20 


o.40!»:|^ 


V1223 Sgr 


0.65 ± 0.04 


0.95 ± 0.05 


1 07+°-"^ 

^■"'-0.09 




1.28(> 0.84) 


RX J2133 


0.91+"" 


0.91 ±0.06 








IGR J17303 


1 06+^-'^ 

i-""-0.14 


1.08 ± 0.07 











This work. iBrunschweiger et alj (l2009l) . ISuleimanov et all (l2005h . iRamsavl ( l200Ch . ICropper et all ( Il999h . lEzuka & Ishidal (Il999h . 



the contrary. iMhlahlo et all (l2007h derived Ki = 74 + 2 km s ' 
from Ha, HB, and Hy lines, a nd used K2 of 360 + 35 km s ' 
reported by IVande Putte et al] (|2003), leading to the WD mass 
= 0.50 + 0.05 Mq which disagrees with the former value 
by more than four times the error range. 

Independently of the radial velocities, the PIN spectrum of 
EX Hya clearly shows a cutoff in the 10-20 keV band (Sect. 
14. Il l, and because of this cutoff together with the Fe line struc- 
ture, the lower WD mass is favored by the fitting. The cutoff 
cannot be explained by a WD of 0.7 Mwd since our model im- 



plies a shock temperature of 30 keV and a similar cutoff-energy, 
and even multiple absorption or reflection cannot fake a lower 
cutoff energy from a spectrum of the high shock temperature. 
The inconsistency with the WD mass of Beue rmann & ReinschI 
(2008) might imply the necessity of some modifications to the 
assumptions of the geometry and the physical processes (Sect. 
12.2b in this object. 
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Fig. 7. Comparison of the WD masses est i mated in the 
present study with those by ("a) ISuleimanov et all (l2005h and (b) 
iBrunschweiger et alJ (l2009h . 



5.2. The WD mass spectrum and its average 



Figure |9] shows the distribution of the WD mass of our sample, 
compared with those of CVs det ermined kinematicall y, mainly 
based on optical/lR spectroscopy dRitter & Kolbl20()3h . The two 
distributions give the mean masses of 0.88 ± 0.25 M© for the 
present study and 0.82 ± 0.23 Mq for the optical sample (errors 
are Icr standard deviations), and the Kolmogorov-Smirnov prob- 
ability is 0.39. Therefore it is inferred with this limited sample 
that the mass spectra of CVs and IPs are statistically not largely 
different. With this distribution, although the systematic uncer- 
tainties (Sect. 15.4b should be examined carefully, we consider 
that the WD masses of IPs are truly distributed, not extracted 
from a single common value with errors caused by the analysis 
method. 

The mean WD mass of our IP sample, 0.88 + 0.25 Mp, can be 
compared with the value of ~ 0.5 Mq inferred bv lKrivonos et alj 
(12007,) based on the analysis of Galactic Ridge X-ray Emission 
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5 Ih 
u 
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TV 



BG 
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M^o (This work) 



Fig. 8. Comparison of the WD masses estimated in the present 
study with those in t he op tical and infra-red wavelengths. 
Reference s are iHellieii (I1997D for XY Ari, IPenningI (11985)) for 
BG CMi, lHemel^ (1 19931) for TV Col. Ilaswell et al. ( 1997) for 
YY Dra, and 'Beuerma nn et alJ 12004") for V1223 Sgr. For EX 
Hya, two different values are presented; 0. 790 ± 0.026 Mq 
and 0.50 + . 05 M p by iBeuermann & ReinschI (l2008h and 
iMhlahlo et al.l (l2007h . respectively (see text). Quoted errors of 
the present result and the previous reports are at 90% confidence 
levels. 



(GRXE) in the hard X-ray band with INTEGRAL/SPl. This com- 
parison is intriguing because IPs are thought to have deep con- 
nection with the origin of the GRXE (especially in the hard 
X-ray band). However, before examining any discussion about 
the discrepancy of the two values, we should consider the in- 
completeness of our sample. The current sample only consists 
of IPs localized to the solar vicinity (< 500 pc) because of 
their intrinsic low luminosity. Because their selection is usu- 
ally flux-limited, it can be biased to IPs with higher WD masses 
which give higher fluxes in the hard X-ray band. In addition, 
the GRXE spectrum may be decomposed into various t ypes of 
X-ray sources over a wide energy range (e.g. Revnivtse v et al.l 
|2006). Among them, RS CVn systems and dwarf novae have 
softer spectra, and if their contribution is considered simultane- 
ously with IPs, the IP component in the GRXE spectrum may 
become harder, leading to a higher WD mass. Then the differ- 
ence of the two values may be reduced. 



5.3. Comparison of Fe abundances 



Fuiimoto & Ishida (1997), lEzuka & Ishi'da' (Il999h . and 
Evans et al.. (2004) estimated metal abundances of IPs by 
analyzing the resolved Fe line complex (6-7 keV) and 
X-ray continua with t he multi-temperature PSR model. 
Fuiimoto & Ishida (Il997h obtained the Fe ab undance of EX 
Hya as 0.6 + 0.2 Zq, which was later revised bv lEzuka & Ishidal 
fll999) to 0.72+Q°2 ^0 toge ther with me asurements of the Fe 
abundances of 12 other IPs. lEvans etaL I (j2004) gave the metal 
abundance of 0.121+°°23 ^0 for Aqr (whilst they did not 
take their result as evidence of genuine sub-solar abundance). 

Since our IP sample is well covered by that of 

lEzuka & Ishidal (Il999l) . we compare our result with theirs in 
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Ta ble 4. Comparison of t he Fe abundances with values reported 
bv lEzuka&Ishidal(ll999h associated with 90% confidence inter- 
vals. 



System 


This work 




Ezuka & Ishida 








ZFe(Zo) 




FO Aqr 


0.10(0.01 - 


0.13) 


0.48(0.33 


- 0.64) 


XY Ari 


0.43(0.40 - 


0.47) 


0.10(0.01 


- 0.23) 


TV Col 


0.49(0.43 - 


0.55) 


0.31(0.22 


-0.41) 


TX Col 


0.47(0.35 - 


0.59) 


0.52(0.30 


-0.78) 


PQ Gem 


0.24(0.19 - 


0.29) 


0.10(0.05 


-0.21) 


EX Hya 


0.56(0.53 - 


0.60) 


0.72(0.60 


- 0.80) 


V2400 Oph 


0.25(0.21 - 


0.27) 


0.39(0.30 


- 0.48) 


AOPsc 


0.44(0.40 - 


0.52) 


0.41(0.30 


-0.55) 


VI 223 Sgr 


0.29(0.27 - 


0.32) 


0.68(0.54 


-0.83) 




0.4 0.6 
Zpj, (This work) 

Fig. 10. Fe abundances determ ined in the present stud y, com- 
pared with the measurements bv lEzuka & Ishidal (Il999h . Source 
names are abbreviated. 



5.4. Limitations and systematic uncertainties of our 
assumptions 

In this section, we discuss the limitations in our assumptions 
of the present model calculation and the WD mass estimation. 
As a whole, these may cause systematic uncertainties, and affect 
our WD mass results (Table|2]i to some extent. Nevertheless, we 
regard the present results as meaningful, because the estimated 
systematic effects are not significantly larger than the statisti- 
cal precision associated with our determinations of Mwd- More 
detailed spectral analyses, which are less subject to the limita- 
tions will be possible with data taken with future high-resolution 
and wide-band instruments, such as the X-ray microcalorime- 
ter and t he hard X-ray imager onboard the ASTRO-H satellite 
(■Takaha shi et al.ll2008 



Table |4] and Fig. [TO] Although the two results consistently in- 
dicate sub-solar abundances, they exhibit no particular correla- 
tion. Considering that the Suzaku data have a much wider spec- 
tral coverage, higher statistics, and a somewhat better energy 
resolution than those from ASCA which Ezuka & Ishida ( 199% 
used, our result s may be seen as more updated. Note also that 
lEzuka & Ishidal (|T999) did not take into account partial covering 
or the reflection component. 

When our 17 objects are divided into seven objects with 
Mwd > 1-0 Mq and the other 10 with Mwd < 1-0 Mq, the 
former and latter sub-samples have averaged Fe abundances of 
0.33 + 0.12 Zq and 0.41 ± 0.15 Zq, respectively. The lack of sig- 
nificant difference ensures that our abundance determinations are 
not strongly affected by the continuum shape. We suppose that 
these different sub-solar values reflect various metal abundances 
of the accreting gas and the companion (non-degenerate) stars. 

Previous studies of summed (unresolved) Fe lines also 
reported sub-so l ar Fe abundances in IPs. Using Ginga, 
iBeardmore et al.l (l2000h reported an underabundance of metals 
by a factor of ~ 2 relative to solar for VI 223 Sgr. This roughly 
agrees with our result, 0.30 - 0.35 Zq. de MartinojtaL (2QQ§i 
gave 0.33+;j{J^ Zq for PQ Gem (RE075 1 -H 14), which is also com- 
parable to ours (0.19 - 0.32 Zq). 



5.4.1 . Dynamics of accreting gas and system geometries 

As detailed in Suleimanov et al. (2005") and'Brunschweig er et alJ 
(2009), the assumption that the accreting gas falls freely from 
an infinite distance is not necessarily warranted, since the gas 
may well form an accretion disk and radiate away part of its 
energy release there. If such a disk is truncated, under plausible 
parameters, at a radius of ~ 3 x 10'^ cm in an IP with a WD 
of 0.9 Mq, the shock temperature, and hence Mwd, would be 
underestimated by ~ 20%. 

Similarly, a fast spin of a WD would also cause the present 
method to underestimate Mwd because of the centrifugal force. 
If corrected for this effect, the WD masses in Table|2]of two rapid 
rotators in our sample, V709 Cas (with a spin period of 312.78 
s) and YY Dra (~ 530 s), could increase by a few tens of percent 
depending on their accretion rates. 

In addition to this, a fast spinning WD might obey a different 
(equatoriaUy elongated) mass-radius r elation (e.g. iNauenber^ 
ll972UGerovannis & Papasotirioul 119971) . However, precise nu- 
merical treatment of this effect exceeds the scope of our present 
study, and we leave it untou ched . 

According to Canalle e t al.l (120051) . the inclusion of the 
dipole curvature in the PSR calculation reduces the calculated 
shock temperature by ~ 10% for Mwd = 1-0 M© under typical 
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Fig. 11. Result of the numerical calculation of the PSR model 
when changing the accretion rate. Shock temperatures are plot- 
ted with solid, long dashed, short dashed, and dotted lines for 
accretion rates of 0.1, 1, 5, and 10 g cm"^ s"'. 



accreting colatitudes, and <2-3% for Mwd ^ 0.8 Mq. At present, 
it is hardly possible to observationally constrain the magnetic 
field structure near the WD magnetic poles. Therefore, we re- 
luctantly neglected the curvature of the WD magnetic field when 
calculating the PSR structure (Sect. l2.2T l. In Table|2] five systems 
have upper WD mass limits below 0.8 M©, and we suppose this 
systematic effect can be safely neglected at least for these light 
systems. 



5.4.2. Accretion rate dependence 

In the spectral fitting, we freely varied the normalization (i.e. ac- 
cretion rate) of the PSR spectral model, which was calculated 
assuming a fixed specific accretion rate of 1.0 g cm"^ s"' (Sect. 
12.3b . Strictly speaking, this means a self-inconsistency, because 
the calculated shock height and the shock temperature should in- 
crease and decrease, respectively, as the accretion rate decreases; 
this, in turn, is because the free-fre e cooUng ra te is proportional 
to the plasma density squared (e.g. lAizulll973h . 

To examine the accretion-rate dependence of our numerical 
result, we performed calculations with several different accre- 
tion rates as shown in Fig.[TT] As noted above, a lower value of 
a gives a lower shock temperature at any Mwd - For higher Mwd, 
difference of the resulting shock temperatures is more or less 
conspicuous. At Mwd = 0.8 M©, we see a difference of a fac- 
tor of 1.13 between the results of a = 0.1 and 1.0 g cm"^ s"', 
and a factor of 1.02 between a = 1.0 and 10 g cm"^ s 
Correspondingly, at Mwd =1-2 Mq, the difference increases to 
factors of 1.33 (a = 0.1 to a = 1.0) and 1.10 (a = 1.0 to a = 10). 
In most systems we analyzed (except for V709 Cas), WD masses 
are below 1.2 Mq, and therefore we consider that this uncer- 
tainty affects the WD mass estimates only slightly, much less 
than ~ 30%. Especially in lower Mwd (^ 0.6 M©), the a de- 
pendence is small, and therefore we presume that the estimated 
shock temperatures and WD masses are considerably secure in 
this mass range. 



5.4.3. PSR base temperature 

In the present study, the PSR temperature was assumed to fall 
to zero at its base although irradiating hard X-rays from the PSR 
may heat the WD photosphere (e.g.^'opper et al. 1998). Indeed, 
non-zero base temperatur es have been s uggested observationally 
by several authors (see. lEzuka & Ishida,1999 , and references 
therein). 

This uncertainty could be avoided by iteratively calculating 
the PSR emission model with a finite base temperature until a 
model fit converges into a self-consistent result. We did not per- 
form this because of the limitation on the computation time and 
the statistical quality of the data. Furthermore, as noted in Sect. 
12.31 emission will become thick near the WD surface, and there- 
fore a more sophisticated spectral model which takes this effect 
into account is needed to be constructed for performing this anal- 
ysis. 



5.4.4. Reflection component 

The spectra of the present IP sample show clear detections of 
the fluorescent Fe line at ~ 6.4 keV. This strongly indicates 
that some of the incident X-rays are emitted from a PSR re- 
flected by cold gas nea r the WD surface or pre-shock region (see 
lEzuka & Ishidalfl999l and references therein). However, in the 
present spectral fitting (Sect. 14.11 ) we did not take into account 
a reflection component from the WD surface explicitly because 
the computation power needed to precisely calculate the reflec- 
tion is unaffordable at this moment. Qualitatively, its effect may 
be represented to some extent by the strongly absorbed compo- 
nent in the partial-covering model, which has a similar spectral 
shape to a reflection component, both contributing mainly in the 

PI N energy band. 

ICropperetalJ(ll998h revealed that an inclusion of a reflec- 
tion component changes the obtained WD masses only slightly 
from those of their partial covering fittings. By analyzing a 
wide-band spectrum of VI 223 Sgr with a reflection component, 
iRevnivtsev et all ( |2004|) derived a WD mass of 0.71 + 0.03 Mq. 
Our result, 0.75 + 0.05 M© well agrees with their value, and thus 
an effect of the reflection component can be considered not as 
very large. A full treatment of a reflection component is left as 
our future work. 



6. Summary 

We numerical ly solved the hydrostatic equations for PSRs of 
mCVs of Sule imanov et al.l C2005), and updated it (Sect. l2!2l i in 
terms of (1) the additional cooling function at higher temper- 
atures and (2) the optically thin thermal emission code, which 
includes Fe emission lines. By constructing the spectral model 
file (Sect. lZ2] i. we fitted 17 IPs observed with Suzaku, and suc- 
cessfully reproduced the Fe line structures in the 6-7 keV band 
and at the same time the thermal cutoff detected in the hard X- 
ray band (Sect. l^Tt . From the best-fit models, we estimated the 
WD masses and the Fe abundances. The average WD mass for 
the present sample were 0.88 + 0.25 Mq. All sources exhibited 
sub-solar Fe abundances, and the average iron abundance was 
O.38+O.14Z0. The WD mass of the newly found IP, IGR J17195- 
4100, was estimated to be 1.03+^22 for the first time. 

The derived WD masses are compared with previous reports 
(Sect. 15.1b . and a rough correlation was confirmed, although re- 
sults on a few systems differ significantly (~ 0.3 - 0.4 Mq). We 
also tentatively compared the mean WD mass with the value in- 
ferred from the analysis of the GRXE in the hard X-ray band 
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(Sect. 15.2b . and found a difference by ~ 0.4 Mq, although this 
may arise from various systematic biases. 
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Appendix A: Spectral fitting result 

In Fig. lA.ll we present the observed spectra overlaid with the 
best-fit model spectra (online version only). 
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